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Abstract 



Results are presented from a search for the pair production of first- and second- 
generation scalar leptoquarks in proton-proton collisions at ^/s = 7 TeV. The data 
sample corresponds to an integrated luminosity of 5.0 fb^^, collected by the CMS de- 
tector at the LHC. The search signatures involve either two charged leptons of the 
same-flavour (electrons or muons) and at least two jets, or a single charged lepton 
(electron or muon), missing transverse energy, and at least two jets. If the branching 
fraction of the leptoquark decay into a charged lepton and a quark is assumed to be 
j6 = 1, leptoquark pair production is excluded at the 95% confidence level for masses 
below 830 GeV and 840 GeV for the first and second generations, respectively. For 
/5 = 0.5, masses below 640 GeV and 650 GeV are excluded. These limits are the most 
stringent to date. 
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1 Introduction 



The structure of the standard model (SM) of particle physics suggests a fundamental relation- 
ship between quarks and leptons. There are many models beyond the SM that predict the ex- 
istence of leptoquarks (LQ), hypothetical particles that carry both baryon number and lepton 
number and couple to both quarks and leptons. Among these scenarios are grand unified the- 
ories null, composite models [3J, extended technicolor models HUH, and superstring-inspired 
models [Zl. LQs are color triplets with fractional electric charge, and can be either scalar or 
vector particles. A leptoquark couples to a lepton and a quark with a coupling strength A, 
and it decays to a charged lepton and a quark with an unknown branching fraction j6 or to 
a neutrino and a quark with branching fraction 1 — /S. In order to satisfy constraints from 
bounds on flavor-changing-neutral-currents and from rare pion and kaon decays ||3l|8l, it is as- 
sumed that LQs couple only to quarks and leptons of a single generation. LQs are classified as 
first-, second-, or third-generation, depending on the generation of leptons they couple to. The 
dominant mechanisms for the production of LQ pairs at the Large Hadron Collider (LHC) are 
gluon-gluon (gg) fusion and quark-antiquark (qq) annihilation, shown in Fig.[lj The dominant 
processes only depend on the strong coupling constant and have been calculated at next-to- 
leading order (NLO) |i9||. The cross section for production via the unknown Yukawa coupling 
A of an LQ to a lepton and a quark is typically smaller. 

This paper reports on a search for pair production of scalar leptoquarks. Several experiments 
have searched for pair-produced scalar LQs but none has obtained evidence for them |IT0l - [l5| . 
This search uses a data sample corresponding to an integrated luminosity of 5.0 fb^^ recorded 
with the Compact Muon Solenoid (CMS) detector during the 2011 proton-proton run of the 
LHC at v^s = 7 TeV. The analysis performed in this paper considers the decay of leptoquark 
pairs into two charged leptons of the same flavor (either electrons or muons) and two quarks; or 
into a charged lepton, a neutrino, and two quarks. As a result, two distinct classes of events are 
selected: one with two high-transverse-momentum (px) electrons or muons and at least two 
high-px jets {lljj) and the other with one high-px electron or muon, large missing transverse 
energy (E^'^®), and at least two high-px jets {Ivjj). 
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Figure 1: Leading order diagrams for the pair production of scalar leptoquarks. 



The CMS detector, described in detail elsewhere l|T6ll , uses a cylindrical coordinate system with 
the z axis along the counterclockwise beam axis. The detector consists of an inner tracking 
system and electromagnetic (ECAL) and hadron (HCAL) calorimeters surrounded by a 3.8 T 
solenoid. The inner tracking system consists of a silicon pixel and strip tracker, providing 
the required granularity and precision for the reconstruction of vertices of charged particles in 
the range < ^ < In in azimuth and \r]\ < 2.5, where the pseudorapidity t] is defined as 
f] = — ln[tan(0/2)], and 9 is the polar angle measured with respect to the z axis. The crystal 
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2 Dataset and object reconstruction 



ECAL and the brass/ scintillator sampling HCAL are used to measure with high resolution the 
energies of photons, electrons, and hadrons for \ f]\ < 3.0. The three muon systems surrounding 
the solenoid cover a region < 2.4 and are composed of drift tubes in the barrel region 
(1^1 < 1.2), of cathode strip chambers in the endcaps (0.9 < < 2.4), and of resistive plate 
chambers in both the barrel region and the endcaps (|//| < 1.6). Events are recorded based 
on a trigger decision using information from either the calorimeter or muon systems. The final 
trigger decision is based on the information from all subsystems, which is passed on to the high- 
level trigger (HLT), consisting of a farm of computers running a version of the reconstruction 
software optimized for fast processing. 

The llj] and Ivj] analyses are performed separately and the results are combined as a function 
of the branching fraction /3 and the leptoquark mass Mlq for first and second generations in- 
dependently. The analysis in all four decay channels searches for leptoquarks in an excess of 
events characteristic of the decay of heavy objects. Various triggers are used to collect events 
depending on the decay channel and the data taking periods as described in Section |2] An 
initial selection isolates events with high-px final-state particles (two or more isolated leptons 
and two or more jets; or one isolated lepton, two or more jets, and large £^'®® indicative of the 
emission of a neutrino). Kinematic variables are then identified to further separate a possible 
leptoquark signal from the expected backgrounds, and optimized thresholds on the values of 
these variables are derived in order to maximize the sensitivity to the possible presence of a 
signal in each decay mode. The variables used in the optimization are the invariant mass of 
jet-lepton pairs (M/y), the scalar sum (Sj) of the pj of each of the final state objects, and either 
the invariant mass of the dilepton pair (M;/) in the lljj channels or E^^^^ in the Ivjj channels. 

Major sources of SM background are Z/ 7*+jets, W+jets processes, and tt. Smaller contributions 
arise from single top production, diboson processes, and QCD multijet processes. The major 
backgrounds are determined either from control samples in data or from Monte Carlo (MC) 
simulated samples normalized to data in selected control regions. 

After final selection, the data are well described by the SM background predictions, and upper 
limits on the LQ cross section are set using a CLg modified frequentist approach [17, 18J. Using 
Poisson statistics, 95% confidence level (CL) upper limits are obtained on the leptoquark pair- 
production cross-section times branching fraction as a function of leptoquark mass (Mlq). This 
is compared with the NLO predictions |9| to determine lower limits on Mlq for = 1 and 
j6 = 0.5. The //// and Ivjj channels are combined to further maximize the exclusion in /3 and 
Mlq, especially for the case /S ~ 0.5, where combining the two charmels increases the sensitivity 
of the search. 

The data and the initial event selection are detailed in Section |2] of this paper, followed by 
a description of the signal modeling and background estimates in Section |3] and Section |4| 
respectively. Section |5] contains the final event selection, and Section |6] describes the systematic 
uncertainties. The results of the search are presented in Section|7]and summarized in Section|8] 

2 Dataset and object reconstruction 

For the first-generation ee;';' analysis, events are required to pass a double-electron trigger or a 
double-photon trigger, with an electron or photon pj > 33 GeV. For the first-generation evjj 
analysis, events are required to pass either a single-electron trigger or a trigger based on the 
requirement of one electron with pj threshold between 17 and 30 GeV, missing transverse en- 
ergy threshold between 15 and 20 GeV, and two jets with pj threshold between 25 and 30 GeV. 
The trigger thresholds vary according to the run period. For the second-generation leptoquark 
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analyses, events are required to pass a single-muon trigger without isolation requirements and 
with a Pt threshold of 40 GeV. For the eejj channel, the trigger efficiency is greater than 99%. 
For the evjj channel, the electron trigger efficiency is measured to be 95%. For the and fivji 
channels, the single muon trigger efficiency is measured to be 92%. 

Electron candidates llT9l are required to have an electromagnetic cluster with pj > 40 GeV and 
pseudorapidity < 2.5 (2.2) for the ee;; {evjj) analysis, excluding the transition region be- 
tween the barrel and the endcap detectors, 1.44 < \r]\ < 1.57. The evjj analysis requires lower 
electron | ?/ 1 to reduce the QCD multijet background, with negligible reduction of the signal ac- 
ceptance. Electron candidates are required to have an electromagnetic cluster in the ECAL that 
is spatially matched to a reconstructed track in the central tracking system in both rj and the az- 
imuthal angle (p, and to have a shower shape consistent with that of an electromagnetic shower. 
Electron candidates are further required to be isolated from additional energy deposits in the 
calorimeter and from reconstructed tracks beyond the matched track in the central tracking 
system. In addition, to reject electrons coming from photon conversion in the tracker material, 
the track associated with the reconstructed electron is required to have hits in all inner tracker 
layers. 

Muons are reconstructed as tracks in the muon system that are matched to the tracks recon- 
structed in the inner tracking system [20|. Muons are required to have pj > 40 GeV, and to 
be reconstructed in the HLT fiducial volume, i.e. with \ri\ < 2.1. In addition, muons must be 
isolated by requiring that the tracker-only relative isolation be less than 0.1. Here, the relative 
isolation is defined as a sum of the transverse momenta of all tracks in the tracker in a cone of 
AR = y(A(^)2 + (A?/)2 = 0.3 around the muon track (excluding the muon track), divided by 
the muon pj. To have a precise measurement of the transverse impact parameter of the muon 
track relative to the beam spot, only muons with tracks containing more than 10 hits in the 
silicon tracker and at least one hit in the pixel detector are considered. To reject muons from 
cosmic rays, the transverse impact parameter with respect to the primary vertex is required to 
be less than 2 mm. 

Jets and are reconstructed using a particle-flow algorithm ||2T]| , which identifies and mea- 
sures stable particles by combining information from all CMS sub-detectors. The £^^® calcula- 
tion uses calorimeter estimates improved by high precision inner tracking information as well 
as corrections based on particle-level information in the event. Jets are reconstructed using the 
anti-A;x [22] algorithm with a distance parameter of R = 0.5. The jet energy is calibrated using 
Pt balance of dijet and 7+jet events pSj . In the ee;; and }i}ijj {evjj and }ivjj) channels, jets are 
required to have pj > 30 (40) GeV, and |?/| < 2.4. Furthermore, jets are required to have a 
spatial separation from electron or muon candidates of AR > 0.3. 

The initial selection of eejj or j-ifijj events requires two electrons or two muons and at least two 
jets satisfying the conditions described above. The two leptons and the two highest-px jets are 
selected as the decay products from a pair of leptoquarks. The invariant mass of the two elec- 
trons (muons) is required to be M// > 60 (50) GeV. To reduce the combinatorial background, 
events with a scalar transverse energy Sj = pj{h) + pjih) + Pt(7i) + Piiji) below 250 GeV 
are rejected. 

For the evjj {jivjj) initial selection, events are required to contain one electron (muon) satisfying 
the conditions described above and at least two jets with px > 40 GeV and E™^^ > 55 GeV. The 
jet Px threshold is higher than that in the dilepton channels to account for jet px thresholds 
in triggers used in the evjj channel. A veto on the presence of extra muons (electrons) is also 
applied. The angle in the transverse plane between the leading px jet and the £^'*® vector 
is required to be A^(£^^^,/i) > 0.5 to reject events with misreconstructed E^^^. In order to 



4 



4 Background estimate 



reduce the contribution from QCD multijet events, the lepton and the are required to be 
separated by A^(£^'®®, /) > 0.8. In addition, events are rejected if the scalar transverse energy 
Sl^' = pj{l)+ E^iss + px(;i) + px(;2) is below 250 GeV. 

The initial selection criteria are summarized in Tabled] 

Table 1: Initial selection criteria in the ee;;, evjj, and }ivjj channels. 



Variable 


eejj 




evji 




Pt(/i) [GeV] 


> 40 


> 40 


> 40 


> 40 


pT(/2)[GeV] 


> 40 


> 40 






\v{h)\ 


< 2.5 


< 2.1 


< 2.2 


< 2.1 


\v{h)\ 


< 2.5 


< 2.1 






Pt(;i) [GeV] 


> 30 


> 30 


> 40 


> 40 


pjih) [GeV] 


> 30 


> 30 


> 40 


> 40 


AR(/,7) 


> 0.3 


> 0.3 


> 0.3 


> 0.3 


£™«« [GeV] 

\A<p{ErMi)\ 






> 55 


> 55 






> 0.5 


> 0.5 


\A(p{Ef^^l)\ 






> 0.8 


> 0.8 


Ml, [GeV] 


>60 


> 50 






Ml^ [GeV] 






> 50 


> 50 


[GeV] 


> 250 


> 250 






Sl^ [GeV] 






> 250 


> 250 



3 Signal and background modeling 

The MC samples for the signal processes are generated for a range of leptoquark mass hypothe- 
ses between 250 and 900 GeV, with a renormalization and factorization scale ]d ~ Mlq. The 
MC generation uses the PYTHIA generator l|2H (version 6.422) and CTEQ6.6 parton distribu- 
tion functions (PDF). The MC samples used to estimate the contribution from SM background 
processes are tt+jets events, generated with MadGraph 125, 26]; single-top events (s, t, and tW 
channels), generated with MadGraph; Z/7*+jets events and W+jets events, generated with 
SHERPA ||27||: VV events, where V either represents a W or a Z boson, generated with pythia; 
QCD muon-enriched multijet events, generated with PYTHIA in bins of transverse momentum 
of the hard-scattering process from 15 GeV to the kinematic limit. The simulation of the CMS 
detector is based on Geant4 Il28l , and includes multiple collisions in a single bunch crossing 
corresponding to the luminosity profile of the LHC during the data taking periods of interest. 

4 Background estimate 

The main processes that can mimic the signature of a leptoquark signal in the II jj channels are 
Z/ 7*+jets, tt, yy+jets, W+jets, and QCD multijets. The Z/ 7*+jets background is determined 
by comparing events from data and MC samples in two different regions: in the region of low 
(L) dilepton invariant mass around the Z boson mass (70 < Mu < 100 GeV for electrons and 
80 < Mil < 100 GeV for muons) and in the region of high (H) mass M// > 100 GeV. The low 
mass scaling factor Rz = Ni/Nf^^ is measured to be 1.27 ± 0.02 for both eejj and }ijijj channels, 
where Ni and Nj^'^ are the number of data and MC events, respectively, in the Z mass window. 
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The number of Z/ 7*+jets events above 100 GeV is then estimated as: 

Nh = Rz<S (1) 

where is the number of MC events with M// > 100 GeV. The estimated number of 
Z/7*+jets events is obtained with the selection criteria optimized for different LQ mass hy- 
potheses and it is used in the limit setting procedure. 

The number and kinematic distributions for the tt events with two leptons of the same flavor is 
estimated from the number of data events that contain one electron and one muon. This type of 
background is expected to produce ee plus and eji final state events with equal probability. 
In the data the number of ee or jiji events is estimated to be: 

trig 

where and ee are the muon and electron reconstruction and identification efficiencies and 
etag ^j.g ^Yie HLT efficiencies to select ee, jiji, and eji events. 

No QCD multijet MC events pass the jijijj final selection. A crosscheck made using a data 
control sample containing same sign muons confirms that the QCD multijet background is 
negligible in this channel. 

For the first-generation LQ analyses the multijet background contribution is estimated from 
a data control sample as follows. The probability that an electron candidate passing loose 
electron requirements additionally passes all electron requirements is measured as a function 
of pj and ?/ in a data sample with one and only one electron candidate, two or more jets and 
low E^'®*. This sample is dominated by QCD multijet events and is similar in terms of jet 
activity to the eejj and evjj analysis samples. A correction for a small contamination of genuine 
electrons passing all electron requirements is derived from MC simulations. The QCD multijet 
background in the final ee;; (ev;;) selection is predicted by applying twice (once) the above 
probability to a sample with two electron candidates (one electron candidate and large E™®*), 
and two or more jets, which satisfy all the requirements of the signal selections. The resulting 
estimate is ~1% (~8%) of the total background for the selections corresponding to the region 
of leptoquark masses where the exclusion limits are placed. 

Contributions to the //;;' background from yV+jets processes and single-top production are 
small and they are estimated using MC simulation. 

In the Ivjj channel, the main backgrounds come from three sources: processes that lead to 
the production of a genuine W bosons such as W+jets, tt, single-top production, diboson pro- 
cesses (WW, WZ); instrumental background, mostly caused by the misidentification of jets as 
leptons in multijet processes, thus creating misidentified electrons or muons and misrecon- 
structed in the final state; and Z boson production, such as Z/7*+jets and ZZ processes. 
The contribution from the principal background, W+jets, is estimated with MC simulation nor- 
malized to data at the initial selection level in the region 50 < Mj < 110 GeV, where Mj is the 
transverse mass calculated from the lepton pj and the Ej"*®. 

In the evjj analysis, the region 50 < Mj < 110 GeV is used to determine both the W+jets and the 
tt normalization factors using two mutually exclusive selections (less than four jets or at least 
four jets with pj > 40 GeV and \i]\ < 2.4) that separately enhance the samples with W+jets and 
with tt events. The results of these two selections are used to form a system of equations: 

Ni = Rtt^i^tt + RwNi,w + Ni,QCD + Ni,o; 

N2 = -Ktt-^fyt + •RwN2,W + N2,QCD + N2,0- 
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5 Event selection optimization 



where N,, N, W/ M,0/ ^(,tt' ^^id Ni^qcd are the number of events in data, W+jets, other MC back- 
grounds, tt-MC, and QCD multijet events obtained from data, passing selection i. The solution 
of the system yields the following normalization factors: Rfj = 0.82 ± 0.04 (stat.)±0.02 (syst.) 
and Rw = 1-21 ± 0.03 (stat.)±0.02 (syst.), where the systematic errors reflect the uncertainties 
on the QCD multijet background estimate. 

In the jivjj analysis the W+jets normalization factor is calculated in the region 50 < Mj < 
llOGeVas: 

N-{N,, + RzNz + No) 
= Nw 

where N, N^-^, Nz, Nq, and Nw are the number of events in data, tt, Z/7*+jets, other back- 
grounds (QCD, diboson, single top production), and W+jets MC samples, normalized to the 
integrated luminosity of the data sample. Rw is measured to be 1.21 ± 0.01 (stat.). To estimate 
the total tt background in the jivjj analysis, the MC prediction is compared with data in a con- 
trol sample enriched with tt events (50 < Mj < 110 GeV, £^'®^ > 55 GeV, and at least four jets). 
The normalization factor of the tt background, 

^ N-{Nw + RzNz + No) 

Ktt = (5) 

is calculated to be 0.88 ± 0.02 (stat.). The calculation of the normalization factors Rw and Rjt is 
repeated using an iterative procedure. 

The contribution from QCD multijet processes after all selection criteria are applied to the }ivjj 
sample is estimated to be negligible. No QCD MC events survive the full selection criteria 
optimized for any leptoquark mass hypothesis. However, as multijet processes are difficult to 
accurately model by MC simulation, several crosschecks are made with data control samples 
to ensure that the QCD multijet background in the jivjj analysis is negligible. The method used 
to determine the QCD multijet background in the evjj analysis is similar to the one used for the 
ee;; channel. 



5 Event selection optimization 

After the initial selection, the sensitivity of the search is optimized by maximizing the Gaussian 
signal significance S/ \/S + B in all channels. Optimized thresholds on the following variables 
are applied for each leptoquark mass hypothesis in the llj] channels: M//, Sj, and My^^. The 
invariant mass of the dilepton pair, M//, is used to remove the majority of the contribution from 
the Z/ 7*+jets background. The variable MJ^'" is defined as the smaller lepton-jet invariant 

mass for the assignment of jets and leptons to LQs which minimizes the LQ - LQ invariant 
mass difference. 

Thresholds on the following variables are optimized for each leptoquark mass hypothesis in the 
Ivjj channels: E^^^, and M;y. A minimum threshold on £^'®* is used, primarily to reduce 
the dominant W+jets background. The variable M/y is defined as the invariant mass of the 
lepton-jet combination which minimizes the LQ - LQ transverse mass difference. In addition, 
a lower threshold is applied on the transverse mass of the lepton (electron or muon) and E^^^ 
in the event, Mj > 120 GeV. 

The resulting optimized thresholds are summarized in Tables |2] and |3] for the 11]] and the Ivjj 
channels, respectively. 
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Table 2: Optimized thresholds for different mass hypothesis of the signal. 



Mlq (GeV) 


250 


350 


400 


450 


500 


550 


600 


650 


750 


850 


900 


> (GeV) 


330 


450 


530 


610 


690 


720 


770 


810 


880 


900 


920 


Mil > (GeV) 


100 


110 


120 


130 


130 


130 


130 


130 


140 


150 


150 


Mfl^ > (GeV) 


60 


160 


200 


250 


300 


340 


370 


400 


470 


500 


520 



Table 3: Optimized thresholds for different mass hypotheses of the Ivjj signal. 



Mlq (GeV) 


250 


350 


400 


450 


500 


550 


600 


650 


750 


850 


S^- > (GeV) 


450 


570 


650 


700 


800 


850 


890 


970 


1000 


1000 


> (GeV) 


100 


120 


120 


140 


160 


160 


180 


180 


220 


240 


M,i > (GeV) 


150 


300 


360 


360 


360 


480 


480 


540 


540 


540 



After the initial selection criteria are applied, the yields in data are found to be consistent with 
SM predictions. Distributions of variables used in the final selection for the eejj, evjj, and 
jivjj analyses are shown in Figs. [2]-[5] 
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Figure 2: ee;;' channel: the distributions of (left) and of M^j for each of the two electron-jet 
pairs (right) for events that pass the initial selection level. The data are indicated by the points, 
and the SM backgrounds are given as cumulative histograms. The expected contiibution from 
a LQ signal with Mlq = 400 GeV is also shown. 
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Figure 3: evjj channel: the distributions of (left) and of Mgj (right) for events that pass the 
initial selection level. The data are indicated by the points, and the SM backgrounds are given 
as cumulative histograms. The expected contribution from a LQ signal with Mlq = 400 GeV is 
also shown. 



The number of events selected in data and estimated backgrounds are then compared at differ- 
ent stages of selection. This information is shown in Tables |4]-[7| for the initial selection and for 
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5 Event selection optimization 
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Figure 4: channel: the distributions of Sj (left) and of M^j for each of the two muon-jet 
pairs (right) for events that pass the initial selection level. The data are indicated by the points, 
and the SM backgrounds are given as cumulative histograms. The expected contribution from 
a LQ signal with Mlq = 400 GeV is also shown. 



> 



CMS 2011 
\/s = 7 TeV 



- Data, 5.0 fb"' 

- W + jets 



'• Other backgrounds 
; M,„ = 400GeV 



jjiij 



1000 1200 1400 1600 1800 2000 
S!f" (GeV) 



> 



d 




1000 1200 1400 1600 
M„:(GeV) 



Figure 5: iivj] channel: the distributions of Sj (left) and of M^j (right) for events that pass the 
initial selection level. The data are indicated by the points, and the SM backgrounds are given 
as cumulative histograms. The expected contribution from a LQ signal with Mlq = 400 GeV is 
also shown. 
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the final selection for each channel separately. 

Table 4: Individual background (BG) sources, expected signal, data and total background event 
yields after the initial (first row) and final selections for the eej) analysis. Other BG includes 
single top,W+jets, 7+jets, and W+jets. Only statistical uncertainties are reported. 



Mlq 


Z+jets 


tt 


QCD 


Other BG 


LQ Signal 


Data 


Total BG 




6234 ± 24 


768 ± 19 


49.59 ± 0.43 


147.6 ± 2.3 




7201 


7199 ± 31 


400 


35.7 ±1.8 


19.1 ±3.1 


0.877 ± 0.022 


3.12 ± 0.56 


487.4 ± 2.2 


55 


58.8 ± 3.6 


450 


15.2 ±1.1 


7.8 ± 2.0 


0.310 ±0.013 


1.92 ± 0.60 


225.6 ±1.0 


26 


25.2 ± 2.3 


500 


6.55 ± 0.70 


2.45 ±1.10 


0.192 ± 0.012 


1.03 ± 0.42 


109.30 ± 0.46 


14 


10.2 ± 1.4 


550 


4.65 ± 0.58 


0.98 ± 0.69 


0.139 ± 0.012 


0.84 ± 0.42 


57.35 ± 0.23 


11 


6.60 ± 0.99 


600 


3.04 ± 0.46 


0.49 ± 0.49 


0.088 ±0.011 


0.72 ± 0.41 


30.95 ± 0.14 


8 


4.34 ± 0.79 


650 


2.14 ±0.38 


0.49 ± 0.49 


0.073 ±0.011 


0.48 ± 0.40 


16.998 ± 0.065 


6 


3.18 ±0.74 


750 


1.04 ± 0.26 


n 000+"-^^ 

U.UUU_gQQ 


0.0092 ± 0.0020 


0.41 ± 0.40 


5.526 ± 0.023 





1 45+0.73 


850 


0.81 ± 0.23 


n 000+"-^^ 

U.UUU_gQQ 


0.00101 ± 0.00022 


0.40 ± 0.40 


1.9679 ±0.0078 





1 91+0.72 
^■^^-0.46 



Table 5: Individual background (BG) sources, expected signal, data, and total background event 
yields after the initial (first row) and final selections for the evj) analysis. Other BG includes 
single top, Z+jets, 7+jets, and VV+jets. Only statistical uncertainties are reported. 



Mlq 


W+jets 


tt 


QCD 


Other 


LQ Signal 


Data 


Total BG 




20108 ± 99 


9301 ± 42 


3267 ±26 


1913 ±53 




34135 


34590 ± 120 


400 


28.7 ±3.6 


17.5 ±1.8 


6.20 ± 0.46 


6.01 ±0.77 


126.01 ± 0.82 


43 


58.4 ±4.1 


450 


19.7 ±2.9 


12.2 ±1.5 


3.01 ±0.31 


4.13 ± 0.44 


68.38 ± 0.43 


29 


39.1 ± 3.3 


500 


13.3 ±2.4 


6.3 ±1.1 


1.72 ± 0.22 


2.80 ± 0.37 


34.70 ± 0.23 


18 


24.2 ± 2.6 


550 


2.98 ± 0.95 


3.38 ± 0.82 


0.65 ±0.10 


1.46 ± 0.26 


16.25 ±0.10 


10 


8.5 ±1.3 


600 


2.45 ± 0.87 


2.33 ± 0.67 


0.57 ±0.10 


1.29 ± 0.25 


9.442 ± 0.056 


6 


6.6 ±1.1 


650 


2.03 ± 0.83 


1.01 ±0.41 


0.335 ± 0.079 


0.76 ± 0.20 


5.202 ± 0.032 


4 


4.14 ±0.95 


750 


1.45 ±0.65 


0.62 ± 0.31 


0.287 ± 0.080 


0.65 ±0.18 


1.851 ± 0.010 


4 


3.01 ± 0.75 


850 


1.22 ±0.61 


0.62 ± 0.31 


0.251 ± 0.078 


0.61 ± 0.19 


0.6973 ± 0.0037 


4 


2.70 ± 0.71 



Data and background predictions after final selection are also shown in Figs. |6]-|9| which com- 
pare St and the best combination for the lepton-jet invariant mass for a signal leptoquark mass 
of 600 GeV in the four decay channels considered. 

6 Systematics Uncertainties 

The uncertainty on the integrated luminosity is taken as 2.2% II29I . 

The statistical uncertainties on the values of Rz (Rw) after initial selection requirements are 
used as an estimate of the uncertainty on the normalization of the Z/7*+jets (W+jets) back- 
grounds. The uncertainty on the shape of the Z/ 7*+jets and W+jets distributions is calculated 
to be 15% (10%) and 20% (11%) for first (second) generation, respectively, by comparing the 
predictions of MadGraph samples produced with factorization or renormalization scales and 
matrix element-parton shower matching threshold varied up and down by a factor of two. 

The uncertainty on the estimate of the tt background in the eejj and fijiji channels is derived 
from the statistical uncertainty of the efiji data sample and the ratio of electron and muon 
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Figure 6: ee;';' channel: the distributions of S'f (left) and of Me; for each of the two electron-jet 
pairs (right) for events that pass the final selection criteria optimized for a signal LQ mass of 
600 GeV. The data are indicated by the points, and the SM backgrounds are given as cumulative 
histograms. The expected contribution from a LQ signal with Mlq = 600 GeV is also shown. 



CMS 2011 
= 7 TeV 



— Data, 5.0 fb"' 

W + jets 

^v$S^ tf 

■rr-r- QCD multijcts 

Other backgrounds 
Mlq = 600 GeV 



2000 2500 
S" (GeV) 



> 



CMS 2011 
\ls = 7 TeV 



- Data, 5.0 tb ' 

- W + jets 
; ti 

QCD multijets 
; Other backgrounds 
Mj^jj = 600 GeV 



lOOO 1200 1400 1600 1800 2000 2200 
(GeV) 



Figure 7: ev]] channel: the distributions of Sj' (left) and of Mgy (right) for events that pass 
the final selection criteria optimized for a signal LQ mass of 600 GeV. The data are indicated 
by the points, and the SM backgrounds are given as cumulative histograms. The expected 
contribution from a LQ signal with Mlq = 600 GeV is also shown. 
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Figure 8: channel: the distributions of Sj (left) and of M^^y for each of the two muon-jet 
pairs (right) for events that pass the final selection criteria optimized for a signal LQ mass of 
600 GeV. The data are indicated by the points, and the SM backgrounds are given as cumulative 
histograms. The expected contribution from a LQ signal with Mlq = 600 GeV is also shown. 



11 



Table 6: Individual background (BG) sources, expected signal, data and total background event 
yields after the initial (first row) and final selections for the jijijj analysis. Other BG includes 
single top, W+jets, and yV+jets. Only statistical uncertainties are reported. 



Mlq Z+jets tt Other BG LQ Signal Data Total BG 



8413 ±46 1028 ±27 197.9 ± 2.8 - 9725 9638 ± 53 



"iOO 45.9 ± 2.4 24.3 ± 4.1 3.54^^;^^ 629.3 ± 4.0 68 73.8+|;^ 

500 10.3 ±1.1 4.2 ±1.7 0.88+°|? 136.70 ±0.86 14 15.3;|J 

550 7.18 ±0.92 2.1 ± 1.2 0.53toil 70.49 ±0.40 9 9.8+J| 

600 4.95 ± 0.74 0.69 ±0.69 0.47l°;f^ 37.39 ±0.22 6 6.1^^;^ 

650 3.76 ±0.64 0.69 ±0.69 0.24+g;f^ 20.56 ±0.13 5 4.7+J;J 

750 2.00 ±0.46 0.00+^;^^ OW^ofs 5.875 ± 0.036 1 2.1+J;^ 

850 1.53 ±0.41 O.OOtoll Omt^fs 2.327 ±0.014 1.6^^;^ 



Table 7: Individual background (BG) sources, expected signal, data, and total background event 
yields after the initial (first row) and final selections for the fivj] analysis. Other BG includes 
single top, Z+jets, and FF+jets. Only statistical uncertainties are reported. 



Mlq 


W+jets 


tt 


Other BG 


LQ Signal 


Data 


Total BG 




23000 ± 160 


13066 ± 61 


2106 ±17 




39287 


38170 ± 170 


400 


30.2 ± 4.2 


25.7 ±2.8 


5.88 ± 0.64 


118.4 ±1.2 


61 


61.8 ±5.1 


500 


17.0 ±3.1 


10.3 ±1.8 


3.40 ± 0.51 


34.02 ± 0.28 


26 


30.6 ± 3.6 


550 


5.6 ±1.7 


4.5 ± 1.2 


2.16 ±0.42 


15.53 ± 0.14 


12 


12.3 ±2.1 


600 


5.4 ±1.7 


3.5 ± 1.1 


1.78 ±0.38 


9.245 ± 0.077 


8 


10.7 ±2.0 


650 


2.8 ±1.1 


1.20 ± 0.62 


1.12 ±0.30 


4.483 ± 0.040 


7 


5.1 ±1.3 


750 


2.8 ±1.1 


0.54 ± 0.38 


0.98 ± 0.28 


1.844 ±0.015 


6 


4.3 ±1.2 


850 


2.7 ±1.1 


0.54 ± 0.38 


0.75 ± 0.23 


0.6934 ± 0.0052 


6 


3.9 ±1.2 
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Figure 9: jivji channel: the distributions of Sj" (left) and of M^y (right) for events that pass 
the final selection criteria optimized for a signal LQ mass of 600 GeV. The data are indicated 
by the points, and the SM backgrounds are given as cumulative histograms. The expected 
contribution from a LQ signal with Mlq = 600 GeV is also shown. 
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7 Results 



reconstruction uncertainties, which is calculated to be 2% and 3%, respectively. A 5.5% (5%) 
uncertainty on the normalization of the estimated tt background in the evjj {jivjj) channel is 
given by the statistical uncertainty on the value of R^i after the initial selection requirements. A 
10% (10%) uncertainty on the shape of the tt background distribution in the evjj {jivj]) channel 
is estimated by comparing the predictions of MadGraph samples produced with factorization 
or renormalization scales and matrix element-parton shower matching thresholds varied up 
and down by a factor of two. 

A systematic uncertainty of 50% (25%) on the QCD multijet background estimate for the eej] 
(evji) channel is estimated from the difference between the number of observed data and the 
background prediction in a QCD multijet-enriched data sample with lower jet multiplicity. 

PDF uncertainties on the theoretical cross section of LQ production and on the final selection 
acceptance have been calculated using the PDF4LHC ||30| prescriptions. Uncertainties on the 
cross section vary from 10 to 30 % for leptoquarks in the mass range of 200-900 GeV, while the 
effect of the PDF uncertainties on signal acceptance varies from 1 to 3%. The PDF uncertainties 
are not considered for background sources with uncertainties determined from data. An uncer- 
tainty on the modeling of pileup interactions in the MC simulation is determined by varying 
the mean of the distribution of pileup interactions by 8%. 

Energy and momentum scale uncertainties are estimated by assigning a 4% uncertainty on 
the jet energy scale, a 1% (3%) uncertainty on the electron energy scale for the barrel (endcap) 
region of ECAL, and a 1% uncertainty on the muon momentum scale. The effect of electron 
energy, muon momentum, and jet energy resolution on expected signal and backgrounds is 
assessed by smearing the electron energy by 1% and 3% in the barrel and endcaps, respectively, 
the muon momentum by 4%, and by varying the jet energy resolution by an ;y-dependent value 
in the range 5-14%. In the Ivjj analyses, the uncertainty on the energy and momentum scales 
and resolutions are propagated to the measurement of £^®®. The effect of these uncertainties is 
calculated for the (minor) background sources for which no data rescaling is applied. For the 
background sources for which data rescaling is applied, residual uncertainties are calculated 
(i.e. relative to the initial selection used to derive the rescaling factor). 

Recent measurements of the muon reconstruction, identification, trigger, and isolation efficien- 
cies using 2. ^ jiji events show very good agreement between data and MC events ||3T1l . A 
~ 1% discrepancy is observed in the data-to-MC comparison of the muon trigger efficiency. 
This discrepancy is taken as a systematic uncertainty per muon, assigned to both signal and 
estimated background. The electron trigger and reconstruction and identification uncertain- 
ties contribute 3% (4%) to the uncertainty in both signal and estimated background for the ee;';' 
(evji) channel. 

The systematic uncertainties, and their effects on signal and background are summarized in 
Table [H] for all channels, corresponding to the final selection optimized for Mlq = 600 GeV. 

7 Results 

The number of observed events in data passing the full selection criteria is consistent with the 
SM background prediction in all decay channels. An upper limit on the leptoquark production 
cross section is therefore set using the CLg modified frequentist approach [17, 18]. A log-normal 
probability function is used to integrate over the systematic uncertainties. Uncertainties of 
statistical nature are described with F distributions with widths determined by the number of 
events simulated in MC samples or observed in data control regions. 
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Table 8: Systematic uncertainties and their effects on signal (S) and background (B) in all chan- 
nels for the MiQ = 600 GeV final selection. All uncertainties are symmetric. 



Systematic 


Magnitude 




eejj 


mil 


evjj 






Uncertainties 


(%) 


S(% 


B{%) 


S(%) 


B(%) 


S(%) B(%) 


S(%) 


B(%) 


Jet Energy Scale 


4 


2 


1 


1 


2 


5 


8.5 


3 


8.5 


Background Modeling 






11 




7 




11 




10 


Electron Energy Scale 


1(3) 


1 


6 






1.5 


4.5 






Muon Momentum Scale 


1 






0.5 


5 






1 


3 


Muon Reco/ID/lso 


1 






2 








1 




Jet Resolution 


(5-14) 


0.5 


0.5 


< 0.5 


0.5 


< 0.5 


2 


< 0.5 


2.5 


Electron Resolution 


1(3) 


0.5 


1 






< 0.5 


1.5 






Muon Resolution 


4 






< 0.5 


5 






< 0.5 


2 


Pileup 


8 


1 


1 


0.5 


< 0.5 


1 


1.5 


1 


1.5 


Integrated Lumir\osity 


2.2 


2.2 




2.2 




2.2 




2.2 




Total 




5 


13 


3 


9 


7 


15 


4 


14 



The 95% CL upper limits on (7 x or cr x 2j6(l — j6) as a function of leptoquark mass are 
shown together with the NLO predictions for the scalar leptoquark pair production cross sec- 



tion in Figs. 10 and 11 The theoretical cross sections are represented for different values of the 
renormalization and factorization scale, }i, varied between half and twice the LQ mass (blue 
shaded region). The PDF uncertainties are taken into account in the theoretical cross section 
values. 

By comparing the observed upper limit with the theoretical cross section values, first- and 
second-generation scalar leptoquarks with masses less than 830 (640) GeV and 840 (620) GeV, 
respectively, are excluded with the assumption that /3 = 1 (0.5). This is to be compared with 
median expected limits of 790 (640) GeV and 780 (610) GeV. 

The observed and expected limits on the branching fraction j6 as a function of leptoquark mass 
can be further improved using the combination of the and Ivjj channels, as shown in Fig- 
ure [12] These combinations lead to the exclusion of first- and second-generation scalar lepto- 
quarks with masses less than 640 and 650 GeV for j5 = 0.5, compared with median expected 
limits of 680 and 660 GeV. 



8 Summary 

In summary, a search for pair production of first- and second-generation scalar leptoquarks has 
been performed in decay channels with either two charged leptons of the same-flavour (elec- 
trons or muons) and at least two jets, or a single charged lepton (electron or muon), missing 
transverse energy, and at least two jets, using 7 TeV proton-proton collisions data correspond- 
ing to an integrated luminosity of 5fb^^. The selection criteria have been optimized for each 
leptoquark signal mass hypothesis. The number of observed candidates for each hypothesis 
agree with the estimated number of background events. The CLg modified frequentist ap- 
proach has been used to set limits on the leptoquark cross section times the branching fraction 
for the decay of a leptoquark pair. At 95% confidence level, the pair-production of first- and 
second-generation leptoquarks is excluded with masses below 830 (640) GeV and 840 (650) GeV 
for fi = 1 (0.5). These are the most stringent limits to date. 
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8 Summary 




Mlq (GeV) (GeV) 

Figure 10: Left (right) frame: the expected and observed upper Hmits at 95% CL on the lepto- 
quark pair production cross section times j6 (2/3(1 — as a function of the first-generation lep- 
toquark mass obtained with the eejj (evji) analysis. The expected limits and uncertainty bands 
represent the median expected limits and the 68% and 95% confidence intervals. The systematic 
uncertainties reported in Table |8] are included in the calculation. The left and middle shaded 
regions are excluded by the current ATLAS limit llSl and CMS limits ||l2l|l3l for j6 = 1(0.5) in 
the eejj (evjj) channel only. The right shaded regions are excluded by these results. The cr^heory 
curves and their bands represent, respectively, the theoretical scalar leptoquark pair production 
cross section and the uncertainties due to the choice of PDF and renormalization/ factorization 
scales. 
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MLQ(GeV) MLQ(GeV) 

Figure 11: Left (right) frame: the expected and observed upper Hmits at 95% CL on the lepto- 
quark pair production cross section times j6 (2/3(1 — j6)) as a function of the second-generation 
leptoquark mass obtained with the (jivji) analysis. The expected limits and uncertainty 
bands represent the median expected limits and the 68% and 95% confidence intervals. The 
systematic uncertainties reported in Table |8] are included in the calculation. The left and mid- 
dle shaded regions in the left frame are excluded by the current ATLAS limit ||T4| and CMS 
limit flT] for /3 = 1 in the mij] channel only. The left shaded region in the right frame is ex- 
cluded by the current ATLAS limit [14J for /3 = 0.5 in the jivjj channel only. The right shaded 
regions are excluded by these results. The CTfj^eory curves and their bands represent, respectively, 
the theoretical scalar leptoquark pair production cross section and the uncertainties due to the 
choice of PDF and renormalization/ factorization scales. 
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Figure 12: Left (right) frame: the expected and observed exclusior\ limits at 95% CL on the 
first- (second-)generation leptoquark hypothesis in the /S versus mass plane using the central 
value of signal cross section for the individual eejj and evjj {jijiji and fivjj) channels and their 
combination. The dark green and light yellow expected limit uncertainty bands represent the 
68% and 95% confidence intervals. Solid lines represent the observed limits in each channel, 
and dashed lines represent the expected limits. The systematic uncertainties reported in Table |8] 
are included in the calculation. The shaded region is excluded by the current ATLAS limits [|14t 
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